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Abstract 

In this paper, we show how to use the analysis of the Lie algebra associated 
with a quantum mechanical system to study its dynamics and facilitate the design 
of controls. We give algorithms to decompose the dynamics and describe their 
application to the control of two coupled spin |'s. 

1 Introduction 

For several quantum mechanical systems subject to a control action, an appropriate model 
is the Schrodinger operator (matrix) equation 

X = -iH(u(t))X, X(0) = 1. (1) 

In this equation, H is called the Hamiltonian and it is an Hermitian matrix function of 
the control u = u(t) and X = X(t) is a unitary matrix for every time t, with 1 denoting 
the identity. Equation (1) is an appropriate model for many quantum phenomena under 
three main assumptions: 1) The quantum mechanical system under consideration can 
be adequately approximated as a system having a finite number of energy levels; 2) The 
interaction with the external environment (decoherence) is negligible; 3) The control, 
which usually represents an appropriately shaped electro-magnetic field, can be treated 
as a classical field (semiclassical approximation). Many physical systems share a model 
of the form (1). Examples are given by systems of particles with spin subject to a control 
magnetic field such as in NMR and EPR, molecular systems where the control is an 
electric field and several implementations of quantum information processing. Here the 
control can be seen as an action allowing us to switch among different Hamiltonians to 
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implement given quantum evolutions (quantum gates). The books [1], [3], [14] present 
physical examples of systems sharing the model (1). 

It is well known (see, e.g, [3], [7]) that the set of operators X reachable for (1) by 
varying the control u is the connected Lie subgroup of U(n) corresponding to the Lie 
algebra £ generated by the set T := {iH{u)\u G U}, that is, the smallest Lie subalgebra 
of u(n) containing T (see the Appendix). The Lie algebra C is called the dynamical 
Lie algebra associated with the system and the associated connected Lie subgroup of 
U(n) will be denoted here by e c . 1 There exists a very simple algorithm to calculate L: 
One starts with a basis of JF, F\, . . . ,F r . If r = n 2 or r = n 2 — 1, one stops because 
C = u(n) or C = su(n), respectively. In this case, e c = U(n) or e c = SU(n) and the 
system is said to be controllable. If this is not the case, one performs the Lie brackets 
of depth 1 [Fj, F k ], j ^ k and select the ones that are linearly independent together with 
{Fi, . . . , F r }, say Di, . . . , D s , if any. Then one performs Lie brackets of depth 2 which 
are Lie brackets of Di, . . . ,D S with the F^s, j = 1, . . . , r and select matrices that are 
linearly independent together with Fi, . . . , F r , D 1 , . . . , D s . One goes on this way until one 
does not find any new linearly independent matrices. The set of matrices thus found is a 
basis of the dynamical Lie algebra C. If the size of this set is n 2 or n 2 — 1, we are in the 
controllable case. Otherwise the system is not controllable. 2 However the Lie algebra C 
gives us information about the nature of the dynamics as we shall see next. 

In section 2 we use the Levi's decomposition of Lie algebras to obtain a decomposition 
of the dynamics for (1). We highlight the simplifications that follow from the fact that 
£ is a subalgebra of u(n). Levi's decomposition is a classical result in Lie algebra theory 
but its impact in quantum control has not been considered before. 3 In section 3, we give 
algorithms to calculate such a decomposition. Algorithms for general Lie algebras are 
known [5] but simplified algorithms can be given in the case of interest here. In some 
cases, our algorithms will be a simplified version of the ones in [5] in some cases different 
algorithms will be given. We shall point out this as we go on. In section 4 we use this 
decomposition for a control problem for two spin \ particles. 

2 Decomposition of Quantum Dynamics 

Every Lie algebra C over the field of reals R is the semidirect sum of a semisimple Lie 
algebra S and the maximal solvable ideal in C, 7Z, called the radical, 4 " that is, 

c = s®n. (2) 

1 Extending this notation, we shall use the notation for the connected Lie group associated with a 
Lie algebra K. See the Appendix for definitions. 

2 See subsection 3.2.1 of [3] for further discussion of this procedure. 

3 One exception is the book [3]. However we shall go beyond what is in this book here by pointing 
out the simplifications in this decomposition in the quantum case and giving explicit algorithms for 
calculation of the decomposition. 

4 Scc definitions in the Appendix. 
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Semidirect sum means that 

[S,S] C S, [S,TZ] C ft, [ft, ft] c ft. (3) 

This is a classical result known as the Levi decomposition (see, e.g., [5]). S is called the 
Levi subalgebra. It is the direct sum of p > 1 simple subalgebras j — 1, . . . ,p, i.e., 

5 = ©5,, (4) 

i=i 

Direct sum means that [5;, <S&] = {0}, when Z ^ 6. In our case, the fact that the dynamical 
Lie algebra £ is a subalgebra of u(n) implies several important simplifications. 

Theorem 1 5 If C C -u(n) then the semidirect sum in (2) is a direct sum, i.e., [S,TZ] = 
{0}, and 1Z is Abelian, i.e., [ft, TZ] = {0}. 

This theorem is a consequence of Lie's theorem (cf., e.g., [6], Corollary A in Section 4.1). 

Theorem 2 (Lie's Theorem) Let TZ^ q< be a solvable Lie algebra of n x n matrices over 
the complex field. Then there exists a change of coordinates (i.e., a similarity transfor- 
mation) to put all the elements in ft (pj in upper triangular form. 

We now give the proof of Theorem 1. 
Proof We first prove that ft is Abelian. We shall consider the field extension (see the 
Appendix and, e.g., [15] for a more in depth discussion) of 1Z to the complex field, ft . 
We shall show that 1Z § is Abelian and this implies 1Z Abelian. In fact ft is solvable 
(Abelian) if and only if TZ is solvable (Abelian). Since TZ^ (rj is solvable, according 
to Lie's Theorem 2 it can be realized as upper triangular matrices. As a consequence, it 
can be written as the sum of two subalgebras: a nilpotent Lie algebra Af, corresponding 
to strictly upper triangular matrices, and an Abelian Lie algebra T, corresponding to 
diagonal matrices in the coordinates indicated in Lie's theorem. Moreover, M is an ideal 
in TZ . Now consider R and P in N ' . Since M is nilpotent, there exists a k > such 
that, ad k R P := [R, [R, [• ■ ■ , [R, P]]\] = where the Lie bracket is taken k times. Now, since 
R is skew- Hermit ian, there is no loss of generality in assuming that R is diagonal, i.e., 
R := diag . . . , i\ n ). Moreover we calculate 

= - Vj, l = 1, n, (5) 

where i := V — 1. From this expression it follows that if ad k x P- x is zero for some k, it 
must be zero for every k and in particular for k — 1. Therefore R and P commute and 
M is Abelian. Consider now N e M and Tel. Since M is an ideal [N, T] is in M and 
therefore [N, [N, T]} = since M is Abelian. The calculation (5) shows that and T 
commute and TZ^ q is the sum of two commuting Abelian subalgebras and it is therefore 
Abelian. 

5 This fact is mentioned in the paper [13] but without a proof. We provide a proof here. 
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The proof that [S,1Z] = use the same calculation (5). Since TZ is an Abelian ideal 
(cf. (3)), ad 2 R P = for every R G TZ and P G S, which, from (5), implies ad R P = 0. □ 

This decomposition of the dynamical Lie algebra C has immediate consequences for 
the Lie group of possible evolutions e c and for the dynamics of the quantum system (1). 
For every control u, the solution of (1) X = X(t) factorizes as 

X(t) = Rf[S j . (6) 

3=1 

Here Sj G e s ' J and R G e 11 and all the factors in (6) commute. Moreover R is itself the 
product of elements belonging to one dimensional subgroups. Write TZ as the sum of one 
dimensional Lie algebras 1Z = 0f =1 7£j, then R = U.f =1 Ri, with Ri G e 7 ^. Controlling 
the system (1) means controlling in parallel the systems Sj = —iHsj(u)Sj, Sj(0) = 
1, j = 1, . . .p, and Ri = -iH m (u)R h Ri(0) = 1, I = 1, . . . , q, where -iH Sj 
and —IHri are the components of —iH(u) in Sj, j = l,...,p, and TZi, I = l,...,q, 
respectively. 

Every, finite dimensional, quantum system (1) has the structure of p + q subsystems 
in parallel of Figure 1. The first p subsystems vary on simple Lie groups for which a 
classification is known [9] , [11]. The remaining q subsystems vary on one dimensional Lie 
groups. The total evolution is the commuting product of the evolutions on the various 
subgroups. To obtain the decomposition of the dynamics, we need to find bases for the 
subalgebras, Sj and TZi, j — 1, . . . ,p, I — 1, . . . q, of C from a basis of C Next, we give 
algorithms for this task. 



3 Algorithms 

The book [5] contains several algorithms for general Lie algebras. In our case, we only 
need to consider are subalgebras of u(n). This allows us in some cases to give new and 
simple algorithms and in other cases to simplify the algorithms of [5]. 

Consider a basis {Li, . . . , L s } of C The calculation of bases of the two subspaces S 
and TZ in (2) is very simple in the case of interest here. It is easily seen that 7Z is the 
center of C, and therefore it is the space of the solutions of the system of s = dim£ 
equations [R, Lj] — 0, j — 1, . . . , s, in the variable R G C. Moreover, using the fact 
that for every semisimple Lie algebra S, S = [S,S], along with [S,7Z] = and that TZ is 
Abelian, we have 

[C,C] = [S®TZ,S®TZ] = [S,S]=S. (7) 

Therefore the set of matrices, [Lj, L k ], j ^ k, spans S. 

The algorithm of [5] to find the solvable radical (cf. Section 2.6 in [5]) first finds the 
product space [C, £} and then the adjoint representations of all elements in a basis of 
C. Then the algorithm solves a linear system of dim[£,£] equations in dim£ unknowns 
obtained by using the Killing form (cf. the Appendix). Then one goes on to calculate the 
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Control u 



5i = -iH sl {u{t))S 1 



S p = -iHsp(u(t))S p 



R! = -iH m (u(t))R 1 



x = nLi Ri Yfj=i Sj 




Figure 1: Structure of a quantum control system. The control u drives simultaneously p+q 
systems on Lie groups (which are simple Lie groups or one-dimensional Lie groups). The 
total evolution X is the commuting product of the evolutions on the various subgroups. 

Levi subalgebra (cf. Section 4.13 in [5]). In our case, we have avoided the calculation and 
storage of the adjoint representation. Moreover, once one knows [C, C] there is no need 
to apply any algorithm to find the Levi's subalgebra S as we have (7). 



The calculation of the simple ideals of S, Sj, j 



,p, is more complicated. We 



follow the path indicated in [5]. A preliminary step is the calculation of the so-called 
primary decomposition of S which is also of interest to understand the structure of S. 

3.1 Calculation of the primary decomposition of S 

The following definition is of interest for general Lie algebras. 

Definition 3.1 A Cartan subalgebra of a Lie algebra S is a nilpotent subalgebra A which 
is equal to its normalizer, that is A = {S G *S| [S*, ^4.] C A}. 

In the case of S C u(n), we have the following 



Proposition 3.2 Every nilpotent subalgebra of u(n) is Abelian. In particular the Cartan 
subalgebra of S is Abelian. 
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Proof. The proof uses the same calculation (5) and argument in Theorem 1 to show that 
if adgX = 0, for some k, and two elements S and X in the Lie algebra, then it must be 
ad s X = 0, i.e., S and X commute. □ 

The following algorithm calculates the Cartan subalgebra for S semisimple and S C 
u(n). 

Algorithm 1 

1. Given the semisimple Lie algebra S, set A = {0}. 

2. Select an element X ^ in S. 

3. Calculate the set of elements in S which commute with X. Call this set V. 

Notice V is also a subalgebra of u(n). & Therefore (just like C above), it has a Levi 
decomposition in its semisimple part which is equal to [D,T>] and the center, C(T>). 
This justifies the next step. 

4. Write V = [V, V] © C{V) where CiV) is the center of V. 

5. Set A = A®C{V). 

6. If [D, V] — Stop and return A as Cartan subalgebra, otherwise set S = [D, V] and 
go to step 2. 

The algorithm converges because at each step S is semisimple and V is a proper 
subspace of S, otherwise S would have an element which commutes with all of S which 
contradicts semisimplicity. To show that the algorithm gives in fact a Cartan subalgebra 
of S we have to show two facts 

1) The resulting A is nilpotent (it is in fact Abelian). 

2) Every S E S which is such that [S, A] C A is an element of A. 

Proof. (Proof of 1) and 2) above) Let us denote by Ak and 2\ the Lie algebras A and 
V obtained after the fc-th step is complete. We show 1) by induction on the steps of the 
algorithm. A\ is definitely Abelian as it is the center of V\. Moreover it commutes with 
T>i. The inductive step shows that these two properties are true at each step. Assume 
they are true at step k — 1. At step k, Ak = Ak-i © C(V k ). However both Ak-i and 
CiT>k) are Abelian and [Ak-\,C(T>k)\ = {0}, because T> k C T>k-i, and we know by the 
inductive assumption that Ak-i commutes with V^-i- This proves that Ak is Abelian. 
Using again Ak = Ak-i © C(V k ), and the fact that V k C V k -i, [Ak^k] = {0} and the 
inductive step is complete. 

To show 2) assume [S, A] C A. Then repeating the argument following (5) [S, A] = 
{0}. Consider X and V 1 at the first step. Write S = Si + S with S G £>i and Si 
in a complement of T>\ in S. Since lei and S commutes with A, then S commutes 

6 The fact that it is a Lie algebra, i.e., [Di,Z?2] G T>, for Z?i,£>2 € T> follows immediately from an 
application of the Jacobi identity. We have [[£>i, D 2 ], X] = -[[D 2 ,X},Di] - [[X,D 1 ],D 2 \ = 0. 
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with X. Since So commutes with X so does Si. So, Si is in V\ and in a complement of 
T>i which implies Si = 0. Therefore S G T>\ and we can write it as S = D\ + A\ with 
Di G [Pi, Pi] and G C(X>i) and therefore G A. At the second step we pick another 
X which turns out to be again in A. Since S commutes with A, D\ is in T> 2 . We can 
then write D 1 = D 2 + A 2 with D 2 G [P 2 , ^2], and A 2 G A Therefore S = D 2 + A 2 + A x . 
Proceeding this way, if the procedure ends after r steps, we have that [D r ,T> r ] = {0} and 
S has the form S = A r + A--1 + • • • + A 1: with A, G A, for every j = 1, r. Therefore 
SgA ' ' □ 



Remark 3.3 We remark that the above algorithm and proof is not derived as a special 
case of the corresponding algorithm in [5] (cf. Section 3.2 in [5]) but it is, to the best 
of the author's knowledge, new. It gives the Cartan subalgebra for the specific situation 
of interest in quantum control. It is simpler than the general algorithm both because it 
involves fewer notions of Lie algebra theory and because it involves fewer operations. 7 

The following definition refers to a general Lie algebra over a general field (cf. [5] 
(Definitions 3.1.1 and 3.1.9)). 

Definition 3.4 A collected primary decomposition of a semisimple Lie algebra 5* with 
respect to a Cartan subalgebra A is a vector space decomposition of the form 

S :=^©Vi©V 2 ©---V r , (8) 

where 



1. The subspaces V/s, j — 1, . . . , r, are invariant under adx, for every X G A, that is 

[AVjjCVj, J = l,...,r. 



2. For every X G A, and every Vj, the minimum polynomial of adx restricted to Vj is 
the power of an irreducible polynomial. 

3. For any two subspaces Vj and 14, there exists an X G A such that the minimum 
polynomials of adx restricted to Vj and Vk are powers of two different irreducible 
polynomials. 

Given A such a decomposition exists and is unique (Theorem 3.1.10 of [5]). 

In our case, for every X, the minimum polynomial of adx restricted to Vj must be of 
the type (A 2 + a 2 ), otherwise adx would have eigenvalues with nonzero real parts and-or 
eigenvalues with geometric multiplicity greater than one. This is not possible because we 
have (in an appropriate basis) ad x = —adx (see Proposition 4.1 in Appendix). 

7 The algorithm of [5] involves, among the other things, the tuning a parameter so that a certain vector 
space has a given dimension. 
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An algorithm to calculate the collected primary decomposition is given below. This 
algorithm was derived by applying the general algorithm presented in [5] to the case 
considered here (cf. Section 4.11 in [5]). 

Algorithm 2 

1. Select an element X G A such that adx has dim C — dim A + 1 different eigenvalues. 
That is, except for the eigenvalue (which has eigenspace equal to A), adx is 
non-degenerate. 

Such elements are called splitting elements and they exist (Corollary 4.11.3 of [5]). 
To find such an X, notice that, if {Ai, . . . , A m } is a basis of A then aaUi , • • • , a^A m all 
commute and they can be simultaneously diagonalized. It is easier then to select real 
coefficients Cj, such that J2JLi CjadA^ has the desired property, and X = Y^JLi c jAj- 
For higher dimensional problems it may be more convenient to use randomized 
algorithms. 

Let X be the selected element, the minimum polynomial is of the form m adx (\) = 
nj=i(A 2 +a|), with the ctj all G R , all different from each other and with one of them 
equal to zero (0 is always an eigenvalue of adx, X being an eigenvector). Moreover, 
from the choice of X being splitting / = dimC ^ dim - A _|_ i anc i m a d x (X) is equal to the 
characteristic polynomial except (possibly) for the power of the monomial associated 
to the eigenvalue 0. 8 

2. Take as Vj the (two-dimensional) eigenspaces associated with the pair of purely 
imaginary eigenvalues corresponding to a 2 . That is Vj = {V G C\(ad x + a 2 l)V = 

n\ • -| dim »S— dim A 

Uj, J 1, . . . , g 

We now prove that the decomposition obtained with the above algorithm is the col- 
lected primary decomposition associated with the Cartan subalgebra A. 
Proof. (Proof of 1, 2 and 3 in definition 3.4). Condition 1 is verified since if X 2 G A, 
and V G Vj, we have (since X and X 2 commute and therefore so do adx and adx 2 ) 
(ad 2 x + a 2 l)adx 2 V = ad X2 {ad 2 x + a 2 l)V = 0. Therefore ad X2 V G Vj as well. Condition 2 
is also verified. Since Vj is two dimensional, the minimum polynomial of adx 2 restricted 
to Vj, for every X 2 G A, must be of the form A 2 + b 2 for real b. Any other form would 
imply that adx 2 has eigenvalues with nonzero real part, which has to be excluded because 
of Proposition 4.1 in the Appendix. Condition 3 is verified taking as element X in A 
precisely the splitting element. 

□ 



8 Notice in particular that the difference between the dimension of C and that of A must be an even 
number. 
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3.2 Calculation of the decomposition in simple ideals 

The primary decomposition is a fundamental tool to explore the structure of a semisimple 
Lie algebra S. Using it, one can directly obtain the decomposition into simple ideals (4). 
The algorithm is given in [5] (cf., Section 4.12) and we report it using our notations. 

Algorithm 3 

1. For every j — 1, . . . , r, calculate the spaces 

oo 

J,:=0a4V„ (9) 

k=0 

where Vj are defined in (8). 

Tj is the smallest ideal containing Vj. ad^Vj, is defined inductively, where adP s Vj = 
Vj, and ad k s Vj = [S, ad^Vj]. 

2. The simple ideals Si, I — 1, . . . ,p in (4) are given by the ideals Tj. Notice that some 
ideals may be coinciding. 

Proof. (Proof of Algorithm 3) The main fact used to justify the algorithm is that the 
primary decomposition (8) is compatible with the decomposition in simple ideals (4). 
This is proved in [5] (Theorem 4.12.1) and it means the following: The Cartan subalgebra 
of S, A, splits in p subalgebras, Aj, j = 1, . . . ,p, i.e., A = @ P j=\Aj, where each Aj is a 
Cartan subalgebra of the corresponding Sj. Moreover, each Vj, j = 1, . . . ,r, in (8) is a 
subspace of one of the Sk, k — 1, . . . ,p. From the latter fact, it follows immediately that 
each Sk, k — 1, . . . ,p, is the smallest ideal generated by one Tj, that is, it is of the form 
(9). In fact, each Tj is contained in some Sk and Sk cannot contain any nontrivial ideal 
other than itself, being simple. □ 



4 Example: Control of two interacting spin ^'s 

Consider the control of two interacting spin \ particles subject to a magnetic field. The 
state of particle 1 (2) lives in a 2-dimensional Hilbert space H\ (H2), so that the state of 
the total system lives in a 4-dimensional Hilbert space Hi <8>7Y 2 - This type of systems are 
of interest for example in quantum computation when one wants to perform operations 
with two quantum bits (cf., e.g., [12]). In the model we shall consider an externally 
applied magnetic field is constant, has nonzero component in the x direction only and 
it only affects the first spin. It is however possible to control the interaction between 
the two spins. There are several ways to experimentally achieve this; see, e.g., [8]. The 
Hamiltonian H{u) in the system's equation (1) has the form 

H = Ui(t)a z <g> a z + u 2 (t)a y <g> a y + a x <g> 1 (10) 
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where <J x ,y,z are the Pauli matrices a x 



which satisfy the commutation relations 

[ia x ,ia y ] = ia z , [ia y ,ia z 



%G X 



[ia z ,ia x 



ia,„ 



(11) 



The solution X of (1) represents the evolution on Hi <g> 7i 2 as the state ip(t) evolves as 
ip(t) = X(t)ip(0) with X(t) the solution of (1). The dynamical Lie algebra is generated 
by {io~ z <E> o~ z ,ia y ® a y , ia x <g> 1}. It is given by 



C := span{ia x <g>l,il <g> a x , ia z <g> a z , ia y 



a y , ia z 



<J y , icr y <S> cr z }, 



(12) 



which is 6-dimensional. As the dimension of the full Lie algebra su(4) is 15, this shows 
that the system is not controllable. The Lie group e c gives the set of reachable evolutions. 
In order to understand the nature of this set and perform control to any possible value in 
it, we apply the analysis developed in this paper. An application of the algorithm for the 
calculation of the (Abelian) radical 1Z discussed at the beginning of Section 3 shows that 
1Z = {0} and that £ is semisimple. To calculate the simple subalgebras, we apply the 
algorithms developed in Section 3. We apply Algorithm 1 to find the Cartan subalgebra 
A. Selecting X = ia x <g> 1 at Step 2 of that algorithm, we find 



A = span{ia x <g> 1, il <g> <r x }. 



(13) 



To find the primary decomposition according to Algorithm 2 we have to select a splitting 
element in A. We write the adjoint representations of ia x <g> 1 and il <g> a x in the ordered 
basis indicated in (12). We have 



adj, 



/0 




















1 









o\ 



1 



-1 




ad, 



il®a x 



(0 














\ 
































1 




















-1 








-1 











V0 








1 





o ) 



Neither ia x 
ad 



51+2jl®cr :I 



\0 -1 0/ 

1 nor il ® o x is a splitting element but ia x <g> 1 + HI ® a x is because 
= adi ax ®i + 2adnt$ ax has dim£ — dirndl +1 = 6 — 2 + 1 = 5 different 
eigenvalues, 0, ±3i, ±i. The eigenspaces corresponding to ±3i and ±1 are given by Vi = 
span{i(a z <g)(Ty + (Ty<g)(j z ),i((T z <g)(T z — (Ty<g)(Ty)}, and V2 = span{i{a z ®a y — a y ®a z ),i{a z ® 
a z + a y ®a y )}. With A in (13) and Vi and V2 above, the primary decomposition is C = A® 
Vi © V2. Using Algorithm 3, we obtain the simple component Lie algebras Sa, Sb, as the 
ideals generated by Vi and V2. We have C = Sa®<Sb, [Sa, Sb] = 0, with the simple Lie 
algebras Sa and Sb given by (cf. (9)) Sa = (Bi° =Q ad^Vi := span{Ai, A 2 , A 3 } , Sb = 
(B%L ad k c V 2 := span{B 1 ,B 2 ,B 3 }, with A 1 := \{a x ®l + l®o x ), A 2 := ^((J z <S>a y +ay<S>a z ), 
A3 := \{a z ®o z -o y ® a y ), and B 1 := \{-a x <g> 1 + 1 ® a x ), B 2 := \{a z ®o y -o y ® a z ), 
B 3 := Ua z ® a z + a y ® a y ). The commutation relations [Ai, A 2 \ = A 3 , [A 2 , A 3 ] = A 1: 
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[A 3 ,A 1 ] = A 2 ([Si,B 2 ] = B 3 , [B 2 ,B 3 ] = B u [B 3 ,B 1 ] = B 2 ), compared with (11), show 
that 1S4 and Sb are both isomorphic to su(2). Writing the Hamiltonian (10) as iH = 
iH^ + iHs, iHj\ = A 3 {u\ — u 2 ) + A 1} iHb = B 3 {u\ + u 2 ) — B 1 , we find that the solution X 
of (1), (10) is the commuting product of the solutions Ua and Ub of the decoupled systems 
tf A = -iHA^ux - u 2 )U A , Ua(0) = 1, and U B = -iH B {ui + u 2 )U B , U B (0) = l,i.e., 
X(t) = UA(t)U B (t) = U B (t)UA(t). We can use u\ — u 2 and u\ + u 2 as independent 
controls to drive Ua and Ub, respectively. Because of the isomorphism between Sa and 
Sb with su(2), both control problems are equivalent to control problems on SU(2) for 
which there exists a large literature. One can for example use the Riemannian symmetric 
space argument of [10] to obtain the minimum (in fact infimum) time control if there is 
no bound on the control. Alternatively one can use an optimal, minimum energy, control 
over a finite time horizon which turns out to be given by elliptic functions as described 
in [4] or a Lie group decomposition technique as in [2] which can be applied when there 
are bounds on the magnitude of the controls. 
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Appendix: Some facts about Lie algebras and Lie groups 

A Lie algebra (see e.g. the textbooks [5], [6], [15]) is a vector space closed under with a 
binary operation called the commutator (x, y) — > [x, y] which is bilinear with respect to 
the vector space sum, skew-symmetric ([x,y] = —[y,x]) and satisfies the Jacobi identity: 
[x, [y, z]] + [y, [z, x]] + [z, [x, y]] = 0. For Lie algebras of matrices that are the ones that 
interest us in this paper, the commutator is taken to be the usual commutator of two 
matrices [A, B] := AB — BA. The Lie algebra u(n) (su(n)) is the Lie algebra over the 
reals of n x n skew- Hermit ian matrices (with trace zero). It has dimension n 2 (n 2 — 1). 
Two Lie algebras C and £ are isomorphic, if there exists a linear one to one and onto 
map <f> : C — > £, such that 4>([A,B]c) = [<j>(A), <f>(B)]c>, for every A and B in £ Here 
[■,-]c and [•, -]c> denote the commutators in C and £ respectively. Given a Lie algebra C 
over the reals, it is possible to define a Lie algebra over the complex field £ which 
is called the field extension of £ q< has the same basis as C and it has the same 
dimension over the complex numbers as C over the reals. The Lie brackets between two 
basis elements give the same result as for the real Lie brackets . Associated with a Lie 
algebra £ of matrices is a Lie group e c which is defined as the set of finite products of the 
form e Ll e L2 ■ ■ ■ e Lf , with L\, . . . , Lf e C, for / > 0. A Lie group is a group in the algebraic 
sense and it is a differentiable manifold. Naturally, the open sets in e c are defined by 
requiring that the map 7T : R n — > e c , 7r(ii, ...,t n ) = e Alh ■ ■ ■ e Antn , for any {A 1: . . . A n } 
basis in C is open, i.e., maps open sets in R ra into open sets in e c . The Lie group 
associated with the Lie algebra u(n) (su(n)) is the Lie group of n x n unitary matrices 
U(n) (n x n unitary matrices with determinant equal to 1, SU(n)). A Lie subgroup 
S of e c is a subgroup in the algebraic which is also a sub-manifold in the sense that the 
topology of S coincides with the one induced by the one of e c . If /C is a subalgebra of 
C, then the Lie group e K is a Lie subgroup of e c provided that the last condition on 
the topology is satisfied. 9 Consider a Lie algebra C, and define inductively the following 

9 This is tacitly assumed anytime we talk about a Lie subgroup in the paper. The result mentioned 
at the beginning of the introduction has to be slightly modified if this last assumption is not verified by 
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sequence of subalgebras C° := C, = [C^ k \C^]. C is called solvable if there exists 

a k such that = {0}. Another sequence is given by £ — £, £fc+i = [AAJ- A 
Lie algebra is nilpotent if, there exists a k such that = {0}. It is called Abelian 
if£i = {0}. If a Lie algebra is Abelian it is nilpotent. If it is nilpotent it is solvable. 
An ideal of £ is a subspace X, such that [£,Z] CI. A Lie algebra C is called simple 
if it has dimension > 1 and it contains no ideals except the trivial ones, {0} and C A 
Lie algebra is semisimple if it is the direct sum (i.e., the sum of vector spaces which 
commute with each other) of simple Lie algebras. Semisimple Lie algebras S have the 
property that [S, S] = S. Consider a subspace M of the Lie algebra C. The normalizer 
of J\f, is the set N C (M) = {L e £\{L,J\f] C J\f}. The center of a Lie algebra C is 
the set C(C) := {L G C\[L, C\ = 0}. It is clear that the center is a subspace of the 
normalizer and both of them (using Jacobi identity) are closed under commutation and 
therefore are subalgebras of C. Let L be a Lie algebra over the field of reals. The adjoint 
representation of the Lie algebra £ is a function ad: C — > -Mdim(£),dim(/:) 10 which maps 
X G C to adx, where adx is a linear map C — > C defined by adxL : = [X, L]. The adjoint 
representation is a representation of the Lie algebra C in that it preserves the basic Lie 
algebra operations. We have ad[x,Y] = [o-dx^ady] and adx is a matrix acting on R dim£ . 
The following fact is used in the paper. 

Proposition 4.1 If £ is semisimple, for every X G C, there exists a basis in C such that 
ad x = —adx- 11 

Proof. Consider the Killing form on C, (Y, Z)k defined as (Y, Z)k '■= Tr(adyadz)- Since 
C is semisimple according to Cartan criterion (see, e.g., [5]) the Killing form is non- 
degenerate. Consider a basis of C, which is orthogonal with respect to the Killing form. 
If we define the transposed of a linear operator O by (Y,O t (Z))k = (OY, (Z))k, it is 
easily seen that the matrix form of the transposed in the orthonormal basis is the usual 
transposed of the matrix form of O. We have, given X, for every Y and Z, (Y, ad x Z)x '■ = 
(ad x Y,Z) K = ([X,Y],Z) K = ([Z,X],Y) K = (Y,[Z,X]) K = -(Y,ad x Z) K . The first 
equality follows from the definition of transposed. Then, we used the cyclic property of the 
Killing form in the third equality and its symmetry in the fourth equality. Summarizing, 
we have (Y,adxZ) k = —(Y,adxZ) k- Since the Killing form is non-degenerate, we must 
have ad x Z = —adxZ and since this is true for every Z the claim follows. □ 



saying that the set of reachable values of X for (1) is dense (in the topology of U(n)) in the Lie group 
associated to the dynamical Lie algebra C. Therefore from a practical point of view there is no difference 
whether e c is or is not a Lie subgroup of U(n). The author wishes to thank Francesco Ticozzi for useful 
discussions on this point. 

w M n ,n denotes the Lie algebra of n x n real matrices with the commutator given by the standard 
matrix commutator. 

n Thc author thanks Richard Ng for working out this proof and for helpful discussions on this paper. 
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